Mechanically induced current and quantum evaporation from Luttinger liquids 
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We investigate transport through a tunnelling junction between an uncorrelated metallic lead and 
a Luttinger liquid when the latter is subjected to a time dependent perturbation. The tunnelling 
current as well as the electron energy distribution function are found to be strongly affected by the 
perturbation due to generation of harmonics in the density oscillations. Using a semiconducting 
lead instead of a metallic one results in electrons being injected into the lead even without applied 
voltage. Some applications to carbon nanotubes are discussed. 
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I. INTRODUCTION 

The influence of static impurities on the transport 
properties of one-dimensional (ID) correlated electron 
systems is well understood [jl], ^] (for a recent review 
on ID physics see e. g. On the contrary, the ef- 

fect of dynamic time-dependent perturbations on trans- 
port from open Luttinger liquids (LLs) has been barely 
addressed in the literature. Yet, investigations into the 
role of dynamic sources in other correlated systems re- 
veal interesting phenomena. So, it was seen in quan- 
tum evaporation (QE) experiments from superfluid he- 
lium Q| that such perturbations can significantly influ- 
ence the structure of the complicated collective particle 
state. A phonon source (which is essentially a dynamic 
perturbation), imbedded into a superfluid is able to ex- 
cite single particles above the energy needed to escape 
from the condensate. In this paper we address the ques- 
tion whether a correlated ID electron system with open 
boundary shows similar physics. Research in this direc- 
tion is particularly interesting because of a possibility to 
make related measurements on single-wall carbon nan- 
otubes (SWNTs). Their electronic degrees of freedom are 
known to be described by the (four-channel) LL model 

An experimental realisation of a periodic perturbation 
can, for example, be achieved by laser radiation which 
causes localised oscillations of the lattice coupled to the 
electron density. Alternatively, it has been realised by 
the oscillation of the nanotube itself || . If the transport 
through the nanotube changes significantly in the pres- 
ence of such perturbation, it could be used to gain infor- 
mation about the cause of the oscillation. This would add 
to already known applications of nanotubes as detector 
or sensor devices p0| . 

The system under investigation is a half-infinite LL 
with an open end. At the boundary it is contacted by 
an electrode which can be metallic, semiconducting or is 
even opened into vacuum (the latter case can be modelled 
by a metallic electrode with the chemical potential sent to 
infinity as in case of the field emission effect O, [l2| ) . We 
are interested in the current-voltage characteristics and 
in the electron energy distribution function (also called 
the total energy distribution function: TED) in the pres- 



ence of a dynamic scatterer inside the LL. The TED is 
proportional to the energy resolved current in the case of 
emission into vacuum. 

In the experimental setups similar to that of Ref. || , 
the nanotube oscillates as a whole. Clearly such oscil- 
lations strongly affect the phase of the electron wave- 
function, but cause little elastic stress and hence only a 
weak back-scattering. It turns out, somewhat surpris- 
ingly, that a dynamic forward scattering perturbation 
non-trivially effects the electron transport. Contrary to 
a back-scattering perturbation, static forward scattering 
contributes neither to the TED nor to the current, so 
that the contribution of the dynamic part can be singled 
out. Because of these reasons, we concentrate here on 
the forward scattering model, which is very simple and 
exactly solvable and yet leads to interesting predictions 
for the electron transport. 

The paper is organised as follows. In the next Section 
we discuss the TED of electrons in the vicinity of the 
open end in the presence of a time-dependent perturba- 
tion. We derive a general formula for the TED. In the 
subsequent Section we derive a similar formula for the 
tunnelling current across the interface to an additional 
electrode. Section [fv| is then devoted to applications of 
the results to different situations and contains discussion 
of the main findings. 



II. TIME-DEPENDENT PERTURBATION IN ID 

Let us begin by recollecting some known results for the 
TED of electrons in the presence of a time-dependent per- 
turbation. Since a static forward scatterer only changes 
the phase of the wave function, it cannot influence 
the TED. As soon as the perturbation becomes time- 
dependent, the scattering process is not elastic any more 
and particles can acquire or lose energy. Hence the TED 
is changed. In the case of a periodic perturbation with 
the frequency O, the electron can lose or gain the energy 
nfl upon scattering, n being an integer. Thereby an infi- 
nite number of equidistant sidebands emerge in the TED 
even in non- interacting systems [|l3| [wj. As a result, 
the TED is distorted in a staircase manner as shown in 
Fig. [|. We stress that the forward scattering is sufficient 
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to cause these effects. 

This can be quantified for our system using the open 
boundary formalism of Ref . Jll| . According to |H| , in an 
open system the electrons can be regarded as chiral par- 
ticles living in an infinite system. Thereby the operator 
ijj(x,t) on the negative half-axis describes the electrons 
moving towards the boundary (right-movers) and for pos- 
itive x it corresponds to electrons moving in the opposite 
direction (left-movers). Since we are dealing with an ex- 
plicitly time-dependent situation the TED n(ui,t) at the 
boundary should be defined using the Wigner represen- 
tation, 

n(w,t)= / dTe lUT {^(0,t-T/2)ij(0,t + T/2)) . (1) 



Because of the open boundary condition imposed on the 
electron wave functions the probability density to find 
any particles exactly at the boundary is, of course, zero. 
However, we assume the TED to be measured on the 
second last site of the underlying lattice model so that 
x in the last formula is actually not but the lattice 
constant ao- 

Since we are dealing with a non-equilibrium situation, 
we resort to the Keldysh formalism [l7|] , in which the 
TED can be conveniently expressed through one of the 
off-diagonal Green's functions (from now on we drop the 
spacial coordinate), 



r g-+(t + T/2 1 t-T/2) 



(2) 



where (we give all Green's functions for future reference), 

g- + (t,t>) = i(^t(f)^(t)), 
g + -(t,t') = -im)^{t'))> 
g—(t,t') = -i(Ty(t)^(t')), 

, 9 ++(M') = -i(fm^ (*'))• (3) 

Here T and T denote the time- and anti-time-ordering op- 
erations, respectively. Notice that these functions are not 
translationally invariant in time because of the explicit 
time dependence of the perturbation U(t). All calcula- 
tions are most transparent in the bosonized representa- 
tion usually used in the theory of LLs . The interact- 
ing electron field ip(t) can then be expressed in terms of 
a free Bose field 4>{t). The physical interpretation being 
that <fi describes the collective low-energy plasmon exci- 
tations in our correlated electron system. At the system's 
boundary we have @: 



J4>(t)/y/g 



\J2-KdQ 



(4) 



This equation contains the Luttinger liquid parameter g, 
which is related to the interaction strength Uo via g — 
(1 + AUq/tt)- 1/2 U|. The mode expansion for the field 
4>(x), in terms of boson creation b q and annihilation b q 



operators, reads 



4>{x) 




e tqx b a 



-a q/2 



(5) 



We assume the system having length L, so that the mo- 
mentum q takes quantised values q = 2im/L, where n is 
an integer. At the end of all calculations we shall send L 
to infinity (any zero-modes omitted in Eq. (^|) will drop 
out in this limit). In this representation, the (perturbed) 
LL Hamiltonian takes the following form, 



<3>0 



+ / dx [A(x, t) + A(-x, t)]p(x) . 



(6) 



where A(x,t) is the forward-scattering time-dependent 
perturbation and v = vp/g is the renormalised sound 
velocity, vf is the Fermi velocity of the non- interacting 
system and p{x) is the particle density operator. Notice 
that we are already working in the chiral representation 
therefore the last integral is performed over the whole 
real axis. For the same reason the scattering potential 
appears in the symmetrised form. In terms of the Bose 
field, the electron density p(x) reads 



p(x) = ^(x)ijj(x) 



1 



d x 4>(x) 



(7) 



The constant contribution (average density) gives rise 
to an overall energy shift and will therefore be dropped 
in what follows. Any additional terms containing back- 
scattering processes are dropped as well. Combining this 
relation with Eq. (|^) we rewrite the perturbation in terms 
of the Bose operators, 



g>0 



u ® = 4 E \/J[ A «W + A_,(t)](&, - b\y 



—aoq/2 



(8) 

where A q (t) is the Fourier transform of the perturbation 
amplitude, A g (t) = J dx e tqx A(x , t) . To evaluate the 
Green's functions all we need is the time dependence of 
the Bose operators b q . In our simple model, the latter 
can easily be obtained as a solution of the corresponding 
equation of motion, which has the following form, 

i^-b\{t) = [blH] = 
= -t, ? 6t (t) _ ^yCT[A g (t) + A_ g (i)]&t e -<W2 . 



The solution of this equation is 
where 



(9) 



/w = \ j T [Aq{T) + A -^)]^ ivq ^ a<,q/2 ■ 



Using this solution we can write down the Keldysh 
Green's function with indices (ij) as products of the bare 
Green's function g$ (t — t') (i. e. without the perturba- 
tion U(t)) and multipliers, which are responsible for the 
breaking of the translational symmetry in time, 

^"(t,0-.9o J (^^)e 4(x(t) - x(t,)) , (10) 
where the phase factors are 

X(t) = I dr [A(-vt, t-r) + A{vt, t-r)]. (11) 

Remarkably, factorisation ( ^p| ) holds for interacting sys- 
tems; it is a consequence of bosonization. 

We are now able to establish the connection between 
the TED in the presence of a dynamic scatterer and that 
of the unperturbed system nn(uj). In order to accomplish 
this we combine Eqs. ([!]), (j|) and ([Toj) . As a result we 
obtain 

dw' 



v^(x(*+r/2)-x(t-T/2)) 



fl2) 



Let us from now on concentrate on a separable har- 
monicperturbation, adequate for describing the setup of 
Rcf. §: 



A (a;, t) = sm[Slt]r](x) . 



(13) 



We define the function T{t) responsible for the geometry 
of the scatterer, 



1 



T(t) = ~ — [v{-vt) + rj(vr)] . 
23 



(14) 



After simple algebra we obtain 



n(w,t)= / ^n (u') J dre^> 
< exp ( — i 2 sin[£lT/2] \z\ sin[Oi + $] 



(15) 



where the constant z and the phase delay <E> between the 
TED oscillations and those of the perturbation are 



/ dr e inT T(t) , 
Jo 

arcsin[Re[z]/|z|] . 



(16) 
(17) 



In realistic systems such as SWNTs the frequency Q is 
expected to be quite high (in the far MHz range) so that 
it is unlikely that the explicit time resolution of the TED 
or other physical properties of the system would be ex- 
perimentally accessible. Therefore it is instructive to in- 
vestigate time-averaged quantities (i. e. averaged over 
the period of the perturbation). We define the averaged 
TED as 



^ )= 2U 



dt n(u>, t) . 



(18) 
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FIG. 1: The TED for a non-interacting system (g = 1) and 
for a weakly interacting LL (g = 0.9) at different oscillator 
strengths. The Fermi energy Ef is set to zero. 



Applying this prescription to Eq. (|l5|) one finds that 



»U'> = / J dre^-^ 

x J [2|2|sin[fiT/2]] 



(19) 



where Jo is the Bessel function of the first kind. The r 
integral can be performed by using the formula 

J dte^ t J {ysm{nt/2}] = 2ir f] S - n) J*(y/2) . 



n— — oo 



As a result one obtains the following expression: 

oo 

n(w)= J2 n (LJ-nQ)J^[\z\}, (20) 

n— — oo 

where z is the only free parameter in this formula. It con- 
tains all essential information about the scatterer prop- 
erties. Hence, the TED of a system with a harmonic per- 
turbation is a weighted superposition of infinite number 
of original unperturbed TEDs centred at energies nfi. 
This result is in accordance with previous findings for 
non-interacting systems [[l3| The important fact is 
that factorisation ( |To| ) continues to hold even for inter- 
acting systems so that formula ( |2fj| ) still applies. 

Due to the presence of the oscillating perturbation the 
actual TED is distorted even for non-interacting systems, 
see Fig. |l|. For a weak scatterer \z\ <C 1 (in the figure 
\z\ is set to 0.5-1.0 in order to make the details of the 
plot clearer) in a non-interacting system, where no(w) 
is just a step function, the particles can barely receive 
energy larger than O upon scattering, so that the TED 
n(u)) only acquires additional step of the width f2 above 
the Fermi energy Ep and a dip of the same width just 
below Ep. At higher \z\ additional sidebands emerge and 
at sufficiently large \z\ the TED is bounded from above 
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only by the actual bandwidth of the host material. The 
latter can be quite high and exceed the work function so 
that particles can leave the system, or to evaporate from 
it in analogy to QE. 

The TED rio(u>) of any unperturbed system is a prod- 
uct of the Fermi distribution function and the energy- 
dependent local density of states (LDOS). Therefore the 
TED n(u>) for an LL differs from that for the non- 
interacting system only by the energy-dependent details 
in vicinity of each step. For example, in Fig. [l] we plotted 
the TED for an LL with interaction parameter g, where 
the LDOS is given by 

^(e) = ^|e|V^. (21) 

C^j, is a constant, which for a non-interacting system co- 
incides with the energy independent LDOS jj[ [l^]. The 
above discussion applies at zero-temperature. It is not 
difficult to generalise the formulas to finite temperatures. 
This is not of immediate interest though as clearly the 
effect of finite temperatures will be to smear out the side- 
bands. We also note that the assumption that the per- 
turbation is harmonic, Eq. (|l]), is not really restrictive. 
The same qualitative behaviour of the TED will persist 
for any periodic perturbation. 

Formula (^o|) has important implications for the field 
emission theory. The TED of emitted electrons above 
the Fermi edge exists only if interactions are present in 
the emitter material and shows a divergent behaviour 
towards Ep jy], [l^]. If, however, we have an additional 
time-dependent perturbation, the high-energy tail can- 
not be accounted for solely by interactions. In such situ- 
ation we expect the field emission (Auger) singularity to 
emerge at every nQ, however, they would be weaker than 
in the case without the time-dependent perturbation. 



III. GENERAL RESULTS FOR THE 
TUNNELLING CURRENT 

We now turn to the calculation of the current-voltage 
characteristics of the system which is contacted by an 
additional electrode at x — 0. The corresponding Hamil- 
tonian, 

H = H LL [4)] + H [c] (22) 
+ 7 [^(0)c(0) + ct (0)^(0)] + U(t) , 

contains a contribution H [c] describing the electronic 
degrees of freedom in the uncorrelated lead (on the right 
of the tunnelling junction). Here 7 is the tunnelling am- 
plitude between the leads. The current operator can be 
derived in the usual way calculating the time derivative of 
the particle number operator on the left or, equivalently, 
on the right side of the junction. Its average should be 
calculated by means of the Keldysh technique because 
of the explicit time dependence of U(t) and the finite 
voltage V which is applied to the junction |l(| |T7j ], 

j(t) = i-y([J (0,t)V(0, t)-^(0,t)c(0, t)]S c )o, (23) 
where Sc denotes the non-equilibrium S'-matrix defined 
on the Keldysh contour C, 

S c = T c exp {-vyj dt' [^(0, t')c(0, t') 

+ct(0,t>(0,t')]}. (24) 

Tc is the contour ordering operation. Performing the sec- 
ond order perturbative expansion in the tunnelling am- 
plitude one obtains the usual expression for the current 



m = 



%r / dt' 



Go "(*' - t)9o "(*,*') - Gr(t - t')g^-{t',t) 



(25) 



+ G+-(t - t)g^ + (t,t') - G^ + {t - t')g+-(t',t) 



We are again using the local Keldysh Green's functions 
defined as in Eq. (pi). The Green's functions on the op- 
posite side of the junction, denoted by Gq (t), are defined 
in a similar way with the change ip ~ * c and depend, con- 



J 



trary to the situation on the left side of the junction, only 
on the time differences. Taking this and the factorisation 
relation ( [l0|) into account we obtain an expression for the 
current similar to Eq. (151), 



j(t) =7 2 / dt'G(t') exp (- i2sin[Qt' /2]\z\sin[nt + $]) , 



the function Q(t) being 

Q{t) = G -M). 9o -(i) - G ~(f)5 ~H) + G+-(-i). 9o -+(t) - Gv+{t)g+-{-t) . 



(26) 



(27) 
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It can be shown, e. g. by doing a perturbative expansion in |z|, that the time-ordered part of this function does not 
contribute to the tunnelling current and can be dropped. 

Again, we expect that it is hardly possible to measure the time evolution of j(t) directly. A more easily accessible 
quantity is the current averaged over one (or more) period of the oscillating perturbation. Performing the average in 
the same way as for the TED we obtain the following general result, 

j = h^ J " M I M<%-{u)9o + («> + an) - G -+(^ + (ln)g+-(w)) . (28) 

n 

We stress again that we have chosen the specific time-dependence of the perturbation, (|]1|), for the sake of clarity. 
A similar formula will hold for any periodic perturbation (with Bessel functions substituted by Fourier coefficients 
relevant for the given perturbation.) 



IV. APPLICATIONS AND CONCLUSIONS 



The junction is biased by a finite voltage V . Setting the chemical potential in the host to zero and the chemical 
potential of the lead to —V, we can immediately write down the Green's functions entering Eq. (p8|), 

G^(lu) = -iQ(u + V)v c , (29) 
Gq+H = ie(-w-V> C) 
9a~{u) = -i©(w)zA/,(w) , 
= i6(-w)zA/,(w) , 

where v c is a constant density of states on the right electrode and v^(u>) is the LDOS in the host which, in general, 
does depend on the energy u in the relevant energy range (set by V and fl). 



Let us first neglect this energy dependence, i. e. take v^{w) 
into (|28] ) we find that because of the sum rule, 



- „(0) 



Then substituting the above Green's functions 



(30) 



the usual linear current-voltage characteristics J19| is restored: 



(31) 



It is remarkable that while the dynamic scatterer profoundly affects the TED, all these contributions completely 
cancel out (for a constant LDOS) in the total current. (A sum rule similar to ( |3C| ) and hence the result ( |3l| ) can 
be shown to hold for arbitrary periodic perturbation.) Thus, as long as the energy dependence of the LDOS can 
be neglected, no trace of the dynamic perturbation can be seen in the current-voltage characteristics. The situation 
changes dramatically when such energy dependence is taken into account. 

It is instructive to start with a model case, when the host metal is uncorrelated but the particle-hole symmetry is 
slightly violated, i. c. the LDOS is weakly energy-dependent in vicinity of the Fermi energy 



i \ (o) i (i) 



(32) 



where is a constant (the measure of the particle-hole symmetry violation). For the current we then obtain 



»<?V 



2v 



(0) 



+ v c v. 



|V/0| 

1 



.2 72 r 



vn £ nJl 

|v/n|+i 



(33) 



Therefore, as soon as the particle-hole symmetry is violated, the current becomes strongly influenced by the source 
of oscillations. Notice that it is essential to have a finite voltage applied across the junction. As is the case for 
the TED, Eq. ((33|) can be thought of as a sum of contributions from different harmonics, that are generated by the 
perturbation U(t). The prefactors J^[|z|] govern the amplitude of n-th harmonic. The argument \z\ of the Bessel 
functions sets a characteristic energy scale and effectively cuts off the sum over n because J„(|z|) as a function of 



6 



i.i 



i 

> 

5 0.9 
-ST 

0.8 



1 2 3 4 5 

v/n 

FIG. 2: Voltage dependence of the differential conductance for a weakly interacting LL (g = 0.9) contacted by a non-interacting 
metallic lead. dj/dV is normalised to g^ 2 v c C^/2-K. For a non-interacting system the plot would be a straight line dj/dV = 1. 

n decreases exponentially for n > \z\. An important consequence is that for \z\ < 1 the perturbation is not strong 
enough to produce higher harmonics and according to Eq. ( p3| ) the linear behaviour of j(V) is restored with an fl 
dependent conductance. 

In real systems, such as SWNTs, the main reason for the energy dependence of the LDOS is, of course, correlations. 
In LLs, the LDOS is strongly energy-dependent, see Eq. (|2l|). Combining Eqs. (p8|), ( p9| ) and ( f2l| ) we obtain the 
expression for the current between the dynamically perturbed LL and a non-interacting lead: 

2 00 

J = 9y v ^{ Jo [M]sgn0O|Vf / 9 + ]T j£[|z|](sgn(V - Sln)\V - fin| 1/s + sgn(V + Qn)\V + iln^ 3 ) } . (34) 

71=1 

Formula (|34]) is the main result of this paper. As expected from the previous discussion, turning off the interactions, 
9=1, restores the linear current- voltage characteristics expected for a tunnelling junction between two uncorrelated 
metallic leads. For the interacting case we predict frequency and amplitude dependent current across the junction. 
Since the SWNTs oscillate at very high frequencies (several 10 10 Hz which corresponds to ~ 10 -4 eV, see e. g. 
Ref. [^)) we would expect that in the transport experiments on oscillating SWNTs even the lowest harmonics should 
be observable. The most interesting features arise in the voltage dependence of the differential conductance dj/dV, 
see Fig. where it is plotted for three different perturbation strengths. At the onset of every harmonic, V/il = 
integer there are pronounced dips at the positions of additional sidebands in the TED. 

Another important setup is when a LL is coupled via a tunnelling junction to a semiconducting electrode (without 
the bias voltage). In this case one of the Green's function is identically zero, Gq + (lu) = as the conductance band 
of the semiconductor (SC) is nearly empty. The other off-diagonal Green's function contains the SC gap A, and a 
constant v s directly connected with the LDOS in the SC in vicinity of the conductance band bottom |20(| , 

Gq~(lu) = -iv s Q(uj - A) . (35) 

Substituting this into Eq. (^) results in 

2 00 

j=gl-u s C^ J2 *|](n«-A)V*. (36) 

n>A/£2 

Contrary to a metallic junction, j is nonzero even in the non-interacting case. As in the previous setups, the current is 
constituted by contributions from different harmonics. However, since the energy of a given harmonic has to overcome 
the SC gap to contribute to the current, the sum starts at n — A/f2 (we assume fi to be positive). The tunnelling 
current shows a very sharp threshold at \z\ = A/fl because for \z\ < A/f2 the contribution of the sum on the right- 
hand-side of Eq. (|3^) is exponentially small. Of course, Eq. (|36|) gives the current immediately after the perturbation 
was switched on since the charging effects across the junction would tend to suppress it. Such process gives rise to a 
finite potential difference across the junction which can be detected. As it is related to the frequency f2 and amplitude 
of the oscillations, such hetero-j unction might be employed as a detector. 




- z=0.5 
-- z=1.0 
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In conclusion, we analysed the influence of a time-dependent forward scattering perturbation on the transport 
properties of correlated one-dimensional electron systems. We derived the exact Green's functions of the system in 
the presence of the dynamic perturbation. This knowledge allowed us to establish the relation between the total 
TED of the electrons in the case with and without the perturbation. In the case of a periodic perturbation with the 
frequency f2 the TED develops sidebands at energies Ep + nfl, n being an integer. The precise energy dependence of 
the sidebands coincides with the behaviour of the local density of states in the particular interacting system without the 
perturbation. We further developed this theory to calculate the current- voltage characteristics of a tunnelling junction 
between an interacting electron system and a metallic or semiconducting electrode and found that the transport is 
strongly affected by the additional dynamic perturbation. 



Acknowledgments 



We would like to thank H. Grabert for a valuable discussion. This work was partly supported by the EPSRC of 
the UK under grants GR/N19359 and GR/R70309. 



[1] C. L. Kane and M. P. A. Fisher, Phys. Rev. B 46, 15233 (1992). 
[2] A. Furusaki and N. Nagaosa, Phys. Rev. B 47, 382 (1993). 

[3] A. O. Gogolin, A. A. Nersesyan, and A. M. Tsvelik, Bosonization and Strongly Correlated Systems (Cambridge University 
Press, 1998). 

[4] M. J. Baird, F. R. Hope, and A. F. G. Wyatt, Nature 304, 325 (1983). 

[5] M. Bockrath, D. H. Cobden, J. Lu, A. G. Rinzler, R. E. Smalley, L. Balents, and P. L. McEuen, Nature (London) 397, 
598 (1999). 

[6] R. Egger and A. O. Gogolin, Phys. Rev. Lett. 79, 5082 (1997). 

[7] C. L. Kane, L. Balents, and M. P. A. Fisher, Phys. Rev. Lett. 79, 5086 (1997). 

[8] C. Dekker, Phys. Today 52, No. 5, 22 (1999). 

[9] P. Poncharal, Z. L. Wang, D. Ugarte, and W. A. de Heer, Science 283, 1513 (1999). 
[10] R. H. Baughman, C. Cui, A. A. Zakhidov, et al., Science 284, 1340 (1999). 
[11] J. W. Gadzuk and E. W. Plummer, Rev. Mod. Phys. 45, 487 (1973). 
[12] A. O. Gogolin and A. Komnik, Phys. Rev. Lett. 87, 256806 (2001). 
[13] J. A. St0vneng and E. H. Hauge, J. Stat. Phys. 57, 841 (1989). 
[14] P. F. Bagwell and R. K. Lake, Phys. Rev. B 46, 15329 (1992). 
[15] M. Fabrizio and A. O. Gogolin, Phys. Rev. B 50, 17 732 (1995). 
[16] L. V. Keldysh, Zh. Eksp. Teor. Fiz. 47, 1515 (1964). 

[17] E. M. Lifshits and L. P. Pitaevskii, Physical Kinetics (Pergamon Press, Oxford, 1981), we are using the notation of this 
book. 

[18] F. D. M. Haldane, J. Phys. C: Solid State Phys. 14, 2585 (1981). 
[19] G. Mahan, Many-particle physics (Plenum press, 1991). 

[20] Strictly speaking one has to set u s (lu) — const lo/\/lu 2 — A 2 , but since this energy dependence does not qualitatively change 
the results, so we drop it. 



